Lattice exciton-polaron problem by quantum Monte Carlo simulations 
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Exciton-polaron formation in one-dimensional lattice models with short- or long-range carrier-phonon inter- 
action is studied by quantum Monte Carlo simulations. Depending on the relative sign of electron and hole- 
phonon coupling, the exciton-polaron size increases or decreases with increasing interaction strength. Quantum 
phonon fluctuations determine the (exciton-) polaron size and yield translation-invariant states at all finite cou- 
plings. 

PACS numbers: 71.35.-y, 02.70.Ss, 63.20.Ls, 71.38.-k 



I. INTRODUCTION 

The binding of electron-hole (E-H) excitations into exci- 
tons (Xs), governing the optical properties of most nonmetal- 
lic materials,' plays a major role in, e.g., organics,^ nanos- 
tructure devices^^ quantum light sources^ Bose-Einstein 
condensation^ and DNA.- 

The coupling of Xs to phonons is widely relevant^iS and 
gives rise to exciton-polaron (X-P) formation, correspond- 
ing to quasiparticles consisting of an E-H pair and a virtual 
phonon cloud. Apart from the essential role of phonons in re- 
laxation processes after optical excitation, lattice-coupling al- 
ters the X radius which determines, e.g., the oscillator strength 
in optics and the overlap of X wave functions required for 
Bose-Einstein condensation. Very recently, a direct obser- 
vation of an exciton-polaron in photoluminiscence spectra of 
quantum dots has been reported.^ 

Examples where X-Ps of intermediate size are clearly 
implicated in current experiments include transition metal 
oxides, such as insulating manganites'" and nickelatesr^ii^ 
though the situation in cuprates is controversial."''^ Another 
important class of materials is conjugated polymers (e.g., 
Ref. In these systems, the well-known approximations 
of small Frenkel or large Wannier-Mott Xs are unjustified, re- 
quiring nonperturbative theories which includes relative E-H 
motion.''* 

Polaron formation is a complex, nonlinear, many-body 
problem which cannot be completely described by renormal- 
ization of effective masses. '^ In particular, the quantum nature 
of phonons — leading to retarded (self-)interaction — has to be 
taken into account. Since polaron physics is governed by lat- 
tice dynamics on the unit-cell scale, the discrete nature of the 
crystal cannot be neglected."" 

The resulting problem of an interacting E-H pair with cou- 
pling to quantum phonons represents a long-standing open 
question in condensed matter physics. Whereas some ex- 
act results are available without phononsr^iS standard meth- 
ods such as perturbation theory, and variational or adiabatic 



approximation; 
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are often of uncertain reliability. 



Furthermore, computational approaches are very demand- 
ing, and we are not aware of any exact results for quantum 
phonons. 

Here we present unbiased numerical results for the quantum 
lattice X-P within a simple E-H model, obtained by means 



of quantum Monte Carlo (QMC) simulations. This method, 
well established in the field of polaron physics, treats all cou- 
plings on the same footing and is not restricted to a specific 
X size or parameter region. Our model study of several dif- 
ferent Hamiltonians yields important results for the effects of 
carrier-phonon interaction on X properties. 



II. MODEL 

Extending previous work^i^l^^ we consider a simple 
model in one dimension (ID) defined by the Hamiltonian 

H = ^ e^j ~ ^ hlhj - ^ Uijhi^eflj.h (1) 
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with long-range Coulomb attraction 
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where Uo > Ui > (i.e., attractive interaction), and long- 
range carrier-phonon interaction 
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Here el (hi) creates an E (H) at site i, and Xi (pi) denotes the 
displacement (momentum) of a harmonic oscillator at site i. 
The fermionic density operators are defined as n^.e — ele^ 
and fii^h — hlh,-. The model parameters are the nearest- 
neighbor E (H) hopping integral (th), the energy of Einstein 
phonons luq (h — 1), the E (H)-phonon couplings (ah), as 
well as the local (extended) Coulomb interaction Uq (Ui). 

We consider a single E-H pair — a situation which 
can be studied experimentally^ — and neglect X cre- 
ation/recombination as well as dynamic screening of the 
Coulomb interaction due to other carriers or lattice polariza- 
tion. Spin degrees of freedom are not taken into account, 
and we assume a tight-binding band structure with s symme- 
try for both E and H, neglecting the existence of a band gap 
(which here only leads to a shift of energies). Of course this 
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model is too simple to make a direct comparison with materi- 
als. Nevertheless, it does describe the physics of a Coulomb- 
bound, itinerant E-H pair whose constituents couple individu- 
ally to quantum phonons and — in the absence of coupling to 
the lattice — captures the familiar crossover from a small to a 
large exciton with increasing bandwidth (see Sec. lIVI) .— 

The exact form of the carrier-phonon coupling is subject 
to X size, screening and material properties.^ We restrict 
our analysis to Holstein- and Frohlich-type interactions well- 
known and understood from polaron physics, and amenable 
to efficient numerical treatment. Important aspects arise from 
the fact that the coupling of E and H to the lattice can either 
be of cooperative or compensating nature. The goal here is to 
obtain a qualitative understanding of the influence of the type 
and range of the lattice coupling, as well as the nonadiabatic- 
ity of the lattice. 

Equation ([T]) allows for different signs of and a^- The 
coefficients fj ^ correspond to a lattice version of the Frohlich 
interaction with longitudinal optical phonons,^'' but yield a 
Holstein coupling to transverse optical phonons for /j ^ = 
6i,j. Since E and H couple to the same phonon mode, we 
consider the symmetric mass case — th — t, and ah — 
aae = <Ja with a — ± and a > 0. We refer to the model 
with local respectively long-range carrier-lattice coupling as 
the Holstein-X model (HXM), respectively, Frohlich-X model 
(FXM). These models capture the interplay of Coulomb at- 
traction, particle motion and coupling to the lattice. 

We introduce the dimensionless parameter A = 
2ep(^^- /|g)/Ty, where sp = a^/2ti;o is the polaron 
binding energy in the atomic limit and W — 4t is the 
bare single-particle bandwidth. The time scales of E/H and 
quantum lattice dynamics are set by the adiabaticity ratio 
7 — ojo/t. The units of energy and length are taken to be Uq 
and the lattice constant, respectively. 

III. METHOD 

The world-line QMC method adapted here can handle long- 
range interactions — notoriously difficult for many other nu- 
merical approaches — higher dimensions, and general fermion 
and phonon dispersion relationsi^i^^ 

From the partition function with discretized inverse tem- 
perature (3 — l/{kBT) and Trotter parameter At = f3/L, the 
fermionic trace can be evaluated using real-space basis states 
{7-^} = {r^, r^}, which define world-line configurations on a 
N X L space-time grid. The path integral over the phonons is 
done analytically, yielding the fermionic partition function 

{'■?} 

P T 

Here carrier-phonon coupling gives the memory function 

F( \ ^ ^0^'^^ V cos[27rT^/£] 

^' AL l-cos[27riy/L] + (woAT)2/2 



with 0(r^ - <0 = Ej y hopping enters via 

1 

Ip{r) = ^ 51 cos(27rfcr/7V)e2^^*''™''(27r/c/7V) 

k=0 

We calculate the X "radius" (see Ref.flTl) 




the kinetic energy 

Ei.in = -t/j2^h+''l^i)^ (8) 

and the binding energies 

E^,^ = Ex{t,Uo,Ui,\)-2E,{t,\) (9) 

and 

^B,A - Exit, Uo, C/i, A) - Exit, Uo, C/i, 0) , (10) 

where Ex (Ee) denotes the X (E) energy. We further study 
the E-H correlation function 

Cehir) = '^ifli^efli+rji) , (H) 
i 

and the E-phonon correlation function 

Cephir) = y^^jfli^eXt+r) ■ (12) 

i 

Computer time ~ (/J/Ar)^ (Ref.lZTl) sets a practical lower 
limit on simulation temperatures. The Trotter error (which 
can be removed by scaling to At = 0, see Ref.l29l) and sta- 
tistical errors limit the accuracy of our QMC results to typi- 
cally 1%, and we use periodic clusters with N = 32. More 
sophisticated QMC approaches to polaron problems, free of 
Trotter errors and finite-size effects,^2!2L^ have been devel- 
oped. Whereas the continuous-time method has recently been 
applied to a similar model,^^ an extension of the diagram- 
matic MC method"^'^"* to the exciton-polaron problem is not 
yet available. Since all three methods are useful only for one 
or two carriers coupled to phonons and hence not applicable to 
more realistic systems, we have chosen the simplest approach 
currently available. 

IV. RESULTS 

To set the stage for the following discussion of lattice ef- 
fects, and to demonstrate that the model defined by Eq. ([T]| 
describes the basic exciton physics, we begin with the case 
A = 0, i.e., no coupling to the lattice. Figure [T] shows exact 
diagonalization (ED) and QMC results for the X size versus 
bandwidth. The zero-temperature ED data for iV = 31 is well 
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FIG. 1 : Exciton radius as a function of bandwidth for A = 
and different values of U\ . Dashed lines correspond to exact ground- 
state results (A'^ = 31). QMC error bars are smaller than the symbols. 



converged with respect to system size. With increasing W , 
there is a crossover from a small, strongly bound Frenkel-X 
with R w (i.e., E and H at the same site) to a larger Wannier- 
Mott-like X with i? > 1. Note that the X is always bound in 
ID.'^ Our parameters do not include the extensively studied 
Wannier-Mott limit, but instead cover experimentally relevant 
intermediate radii. The crossover point (VF/i7o ~ 2) sep- 
arates regions with opposite dependence of i? on t/i.— The 
QMC results are overall in good agreement with T = ED 
data, with finite-temperature effects being most noticeable for 
t/i = 0. 

In the sequel, we restrict ourselves to the wide-band case 
WjU^ = 3.2, highlighted in Fig.[I] for which R(\ = 0) > 1. 
As this work is concerned with phonon effects, we focus on 
the dependence on A and 7, and only consider Ui/Uq — 0.75 
and /3Uo = 15. 

Discussing carrier-quantum-phonon interaction, it is crucial 
to distinguish between a — + and — , as well as between the 
adiabatic (slow lattice, 7 ^ 1) and the non-adiabatic (fast 
lattice, 7 1) regime, taking 7 — 0.4 respectively j — 4. 
We begin with the HXM in the adiabatic regime and <t = +. 

Figure |2ja) shows i? as a function of A. With increasing 
coupling, there is a gradual crossover to a small X-P due to 
the increasingly strong phonon-mediated attractive interaction 
between E and H. The E- and H-polarons tend to maximize 
both the Coulomb and the lattice energy by forming a state 
with small R, but compete with the kinetic energy of the sys- 
tem which decreases with increasing A (Fig. O. Similar to 
the bipolaron problem with U = Q,E- and H-polarons form a 
(phonon) bound state at any A > in ID. Most notably, there 
is no discontinuity at a critical A, a common misconception 
due to earlier variational treatments, as quantum lattice fluc- 
tuations give rise to a translational invariant Bloch-like X-P 
state. 

The crossover is also reflected in a reduced X mobility, and 
in a more negative X binding energy [Fig. |5a)]. With the 
present method, dynamic quantities such as the effective exci- 
ton mass cannot be accurately calculated. An alternative ob- 
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FIG. 2: (Color online) QMC results for _R as a function of A for (a) 
the HXM (see text) and (b) the FXM for both the adiabatic (7 = 0.4) 
and the nonadiabatic (7 = 4) regime, as well as cr = ±. Here and 
in subsequent figures /3t/o = 15, W IU<i = 3.2, C/i/C/o = 0.75, and 
lines connecting data points are guides to the eye. 



servable which to some degree (see, e.g., Ref. 35) measures 
the mobility is the kinetic energy shown in Fig. [3] In addition, 
the E-H and E-phonon correlation functions in Fig. [Sja), al- 
ways positive for ct = +, fall off quickly with r, indicating 
that the X-P is a quasiparticle consisting of a tightly-bound 
E-H pair with a strongly localized suiTounding lattice distor- 
tion. Such a state is similar to the Frenkel limit considered in 
Ref.lM 

The nonadiabatic regime 7^1 mainly differs by a weaker 
dependence on A (the important coupling parameter is ep /wq, 
see below). The results for the FXM exhibit qualitatively 
the same tendencies, but the long-range interaction generally 
leads to larger E- and H-polarons and a larger X-P. 

Turning to the case cr = — , we briefly discuss the different 
polaron ground states in the Holstein and the Frohlich model. 
The Holstein model in ID exhibits a crossover from a large 
polaron to a small polaron (with a predominantly onsite lat- 
tice distortion) with increasing A. For 7^1, the latter occurs 
near A « 1, whereas for 7^1 the condition is epjujQ > 1 
(A = 2 for 7 = 4)."^^ In contrast, the Frohlich polaron re- 
mains large (spatially extended lattice polarization) even for 
strong coupling. While for a — + the bipolaron effect dom- 
inates, these differences have a major impact for a — ~ where 
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FIG. 3: (Color online) Kinetic energy iSkir 




FIG. 4: (Color online) Binding energies Eb,u and Eb.\ (inset). 



the Coulomb-bound E- and H-polarons remain separated with 

In Fig. 121 strikingly different to a = +, R initially increases 
with increasing A, i.e., the X-P is larger for stronger coupling. 
In the HXM [Fig.|2a)], R takes on a maximum at A « 0.7 and 
approaches R = 1 for A > 1, whereas in the FXM R increases 
monotonically and saturates at large R in the strong-coupling 
regime [Fig.|2jb)]. Accordingly, the kinetic energy in Fig.[3]is 
much larger compared to cr = +, but is eventually reduced for 
large A in the HXM. The binding energy E'b.c/ — > with in- 




FIG. 5: (Color online) Electron-hole [Ceh{r)] and electron-phonon 
[Ceph{r), inset] correlation functions for A = 1. 



creasing coupling in both models (Fig.|4|i, whereas i?B,A (see 
insets) — ^related to X-P effects — remains clearly negative. 

The (initial) increase of the radius with increasing A for the 
(HXM) FXM (Fig. ^ is due to the fact that the X-P loses 
lattice energy if the compensating displacement clouds sur- 
rounding E and H overlap. Therefore, the E- and H-polarons 
optimize R to achieve maximum Coulomb energy and mini- 
mum phonon-cloud overlap. The resulting average distance R 
depends on the size of the individual (E and H) polarons. For 
the HXM, polarons are large for A < 1, leading to large values 
of R in Fig. 12 a), and become small for A > 1, causing the de- 
crease of i? ^ 1 in the strong-coupling regime. In the FXM, 
polarons remain large for all A, leading to large values of R 
even for strong coupling [Fig. |2b)]. The larger radius in the 
non-adiabatic HXM in Fig. |2a) as compared to 7 ^ 1 is due 
to the much larger polaron kinetic energy.^^ A discontinuous 
dissociation of the X-P with increasing A has been discussed 
in a continuum model with acoustic phonons and a = — 

From the cr = — results for the E-H correlation function in 
Fig. |5]we see that the E-H separation is small in the HXM, 
whereas the pair is spread out in the FXM. For the HXM with 
7 = 0.4, we find a charge-transfer X-P with E and H mainly 
on neighboring sites. Note that for tr = — , Coulomb and 
carrier-phonon interaction have swapped roles as compared 
to bipolaron formation where the lattice-coupling creates an 
attractive interaction that competes with Coulomb repulsion.— 
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Turning to the E-phonon correlation functions in Fig. |5] we 
have Ceph > for small r, but C^ph < at larger distances as 
a result of the opposite distortions created by the hole. Again, 
the extent of the distortions is much larger for the FXM. 

Real materials will require more detailed modeling, but we 
note that charge-transfer Xs in oxides will be better modeled 
hy a — — (for breathing modes), whereas the characteristic 
case for a direct X in a neutral semiconductor would be ct = 
+. 



V. CONCLUSIONS 

In summary, we have studied the exciton-polaron problem 
with quantum phonons by Monte Carlo simulations. Our sim- 
ple models encompass short- and long-range carrier-phonon 
interaction of either the same or opposite sign for electron and 
hole. There are no sharp transitions with increasing carrier- 
phonon coupling, and for couplings of opposite sign the ex- 
citon radius increases with increasing coupling as a result of 
polaron-polaron repulsion. To capture this effect (depending 
on polaron size which is affected by nonadiabaticity) relative 



electron-hole motion and quantum phonon fluctuations must 
be taken into account. Our findings are expected to be impor- 
tant in materials with relatively small excitons such as organ- 
ics and transition metal oxides, although more realistic models 
will have to be studied for direct comparison. 

The present study motivates future work in a number of dif- 
ferent directions, including more general models with respect 
to band structure, phonon dispersion, spin dependence, disor- 
der, or dimensionality, and many-X-P as well as X-polariton 
problems. To this end, the development of more elaborate 
numerical approaches is highly desirable, permitting investi- 
gations of spectral properties routinely studied experimentally 
or even time-resolved studies of X formation.ii^ 
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